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Remote  Sensing  of  Atmospheric  Crosswinds  by  Utilizing 
Speckle-Turbulence  Interaction  &  Optical  Heterodyne  Detection 


Summar\ 


Speckle-turbulence  interaction  can  be  utilized  to  measure  the  vector  wind 
in  a  plane  perpendicular  to  the  line  of  sight  from  a  laser  transmitter  to  a 
target.  A  continuous  wave  source  of  around  one  watt  and  operating  at  10.6 
microns,  in  conjunction  with  an  optical  heterodyne  receiver,  has  been  used  to 
measure  atmospheric  winds  along  horizontal  paths.  A  theoretical  basis,  the 
experimental  apparatus,  processing  techniques,  and  experimental  results  are 
presented.  The  technique  has  been  demonstrated  for  remote  sensing  of 
atmospheric  winds  along  horizontal  paths;  but  also  has  potential  for  global 
remote  sensing  of  atmospheric  winds  and  for  on-board  wind  shear  detection 
systems  for  aircraft. 

Research  Progress 

o  Designed,  constructed  and  tested  a  10.6  microns,  cw,  optical  heterodyne 
transmitter  receiver  system  for  measuring  atmospheric  winds 

o  Experimentally  confirmed  that  speckle-turbulence  interaction  can  be  used 
to  accurately  measure  atmospheric  winds 


Made  measurements 
turbulence  levels 
good  results 


at  target  ranges  of  500  m,  1000  m  and  1300  m  and  at 
from  10'^^  m"^^^  to  as  low  as  7.4  x  10"^®  with 


0  Discovered  a  novel  method  of  isolating  the  optical  local  oscillator  from 
the  optical  transmitter  in  optical  heterodyne  systems 

o  Developed  three  new  methods  of  processing  the  signals  to  obtain  the  wind 

Additional  detail  is  contained  in  the  three  progress  reports  for  the 
subject  contract,  the  publications  and  reports  listed  later  in  this  final 
report,  and  in  Appendix  A  which  contains  a  preprint  of  a  technical  paper  on 
this  work. 
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ABSTRACT 

Speckle-turbulence  interaction  can  be  utilized  to  measure  the  vector 
wind  in  a  plane  perpendicular  to  the  line  of  sight  from  a  laser  transmitter 
to  a  target.  A  continuous  wave  source  of  around  one  watt  and  operating 
at  10.6  microns,  in  conjunction  with  an  optical  heterodyne  receiver,  has 
been  used  to  measure  atmospheric  winds  along  horizontal  paths.  A 
theoretical  basis,  the  experimental  apparatus,  processing  techniques,  and 
experimental  results  are  presented.  The  technique  has  been  demon- 


Introduction 


Speckle-turbulence  interaction  has  the  potential  for  allowing  single  ended  remote 
sensing  of  the  path  averaged  vector  crosswind  in  a  plane  perpendicular  to  the  line  of 
sight  to  a  target.  If  a  laser  transmitter  is  used  to  illuminate  a  target,  the  resultant 
speckle  field  generated  by  the  target  is  randomly  perturbed  by  the  atmospheric  tur¬ 
bulence  as  it  propagates  back  to  the  location  of  the  transmitter-receiver.  When  a 
crosswind  is  present,  this  scintillation  pattern  will  move  with  time  across  the  receiver 
aperture;  and  consequently,  the  time  delayed  statistics  of  the  speckle  field  at  the 
receiver  are  dependent  on  the  crosswind  velocity. 

A  continuous  wave  (cw)  laser  transmitter  of  modest  power  level  in  conjunction 
with  optical  heterodyne  detection  can  be  used  to  exploit  the  speckle-turbulence  interac¬ 
tion  and  measure  atmospheric  winds.  The  use  of  a  cw  transmitter  at  10.6  microns  and 
optical  heterodyne  detection  has  many  advantages  including  the  availability  of  compact, 
reliable  and  inexpensive  transmitters;  better  penetration  of  smoke,  dust  and  fog;  stable 
output  power;  low  beam  pointing  jitter;  and  relatively  simple  receiver  electronics.  In 
addition,  with  a  cw  transmitter,  options  exist  for  processing  the  received  signals  for  the 
crosswind  that  do  not  require  a  knowledge  of  the  strength  of  turbulence. 

It  should  be  emphasized  that  the  system  to  be  described  does  not  operate  on  the 
same  principles  as  a  doppler  lidar  remote  wind  sensing  system  [1].  The  doppler  systems 
use  the  aerosols  in  the  atmosphere  to  scatter  some  of  the  transmitted  energy  from  a 
coherent  pulsed  laser  system  back  to  a  receiver  where  the  doppler  shift  is  used  to  meas¬ 
ure  the  wind  magnitude  along  the  line  of  sight.  The  speckle-turbulence  system  uses  .1 
hard  target  as  the  scattering  medium  and  measures  the  path  averaged  vector  wind  in  a 
plane  perpendicular  to  the  line  of  sight.  The  two  approaches  compliment  each  other 
and  each  has  certain  advantages.  A  significant  advantage  of  the  doppler  system  is  that 


path  resolved  wind  is  easily  obtained.  There  are  potential  methods  such  as  crossed 
beams  and  multiple  detectors  for  achieving  some  path  resolution  with  a  speckle- 
turbulence  system;  however,  it  would  greatly  complicate  the  system.  The  advantages  of 
a  speckle-turbulence  system  are  the  ability  to  measure  the  vector  wind  in  a  plane  per¬ 
pendicular  to  the  line  of  sight  and  very  simple  equipment  and  data  processing.  By 
using  optical  heterodyne  detection,  only  a  watt  or  two  of  optical  power  is  needed;  and 
consequently,  laser  diodes  could  be  used  as  the  transmitter  source.  This  should  have  a 
positive  impact  on  portability  and  reliability. 

An  example  of  how  well  the  speckle-turbulence  system  measures  the  wind  is  shown 
in  Figure  1.  The  target  was  at  1  km  and  the  turbulence  level  was  intermediate  at 
2.37  X  10“^"*  Forty  measurements  with  two  and  one-half  second  averages  were 

made  and  the  data  was  processed  using  a  new  method  called  the  Binary-Z  Log  Ratio 
method  that  will  be  described  later.  There  is  some  fluctuation  about  the  mean  due  to 
the  speckle;  but  the  RMS  error  is  only  0.33  m/s.  The  mean  value  of  2.01  meters/second 
averaged  over  the  forty  measurements  compares  well  with  the  in  situ  measurement  of 
2.09  meters/second.  Additional  data  will  be  shown  in  a  later  section  including  some 
data  that  demonstrates  good  operation  at  turbulence  levels  as  low  as  7.4  x  10~'" 

Background 

Ajialytic  formulations  for  the  first  and  second  order  statistics  of  the  received  inten¬ 
sity  have  been  developed  for  the  case  where  a  diffuse  target  is  illuminated  by  a  TEMqo 
monochromatic  laser  beam  [2,3].  They  properly  account  for  the  effect  of  atmo-^phr-ri' 
turbulence  both  for  propagation  of  the  laser  beam  from  the  transmitter  to  the  target 
and  propagation  of  the  target  induced  speckle  field  back  to  the  receiver.  The  exact  for¬ 
mulations  are  quite  complex  and  difficult  to  interpret.  However,  for  longer  wavelength 


sources  such  as  10.6  microns,  simple  formulations  based  on  the  joint  Gaussian  assump¬ 
tion  are  valid  for  turbulence  levels  corresponding  to  Rytov  variances  less  than  about 
0.02  and  are  much  easier  to  deal  with  initially. 

From  previous  work,  the  time  lagged  covariance  (TLC)  function  of  the  intensity 
using  the  joint  Gaussian  assumption  is  given  for  the  focused  case  by  [2] 

Ci(P,Vt)  =  <  (l(P2,t2)  -  <I>)  (l(P„ti)  -<!>)> 

_p2  * 

=  <I>^  exp  - r - 5.82  Lk^  J  Cn(w)l(l— w)P— dw  (1) 

2ao  0 

where  Pj  designates  the  location  of  a  detector  in  the  receiver  plane, 

P  =  P2  -  Pi 

T  =  t2  -  tj 

tto  =  Transmitter  Beam  Radius 

Cn  =  Strength  of  Turbulence 
L  =  Path  Length 

V  =  Vector  Wind  Velocity 
k  =  Wave  Number 
and 

w  =  Normalized  Path  Length  from  the  Transmitter-Receiver  to  the  Target 

As  can  be  seen  from  Eq.(l),  if  two  detectors  are  separated  by  a  vector  distance  P.  fin 
TLC  can  be  measured;  and  the  measured  quantity  will  be  a  function  of  the  vector  wind 
velocity.  In  the  next  section  it  will  be  seen  how  to  exploit  this  relationship  and  measure 
atmospheric  winds. 
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Processing  Techniques 

There  are  several  well  known  processing  techniques  that  should  be  useful  in  deter¬ 
mining  the  wind  from  the  received  signal.  Some  of  these  techniques  are  illustrated  in 
Figure  2.  The  Briggs  method  measures  the  time  delay,  Tq,  at  which  the  autocovariance 


and  TLC  curves  cross.  The  delay  to  peak  method  measures  the  time  delay,  Tp,  where 


the  TLC  curve  peaks.  The  width  of  the  autocovariance  method  measures  the  time 
delay,  at  which  the  autocovariance  curve  decreases  to  67%  of  its  peak  value.  The 
slope  method  measures  the  slope  of  the  TLC  at  zero  time  delay.  In  order  to  process  the 
signals  for  the  wind,  the  relationships  between  the  time  delays,  slope,  and  the  wind 
must  be  derived  using  Eq.(l).  It  w’ill  be  found  that  the  delay  to  peak  and  Briggs  tech¬ 
niques  yield  the  wind  directly;  whereas  the  width  of  the  autocovariance  and  slope  tech¬ 
niques  either  require  a  knowledge  of  p^,,  the  transverse  phase  coherence  length,  or  must 
be  used  in  conjunction  with  one  of  the  other  techniques  to  determine  the  wind  and  p,,. 

The  time  delay  to  the  peak  of  the  TLC  function  can  be  found  from  Eq.(l)  by 
differentiating  with  respect  to  t  and  setting  the  resulting  equation  to  zero.  Before 
proceeding,  Eq.(l)  will  be  examined  in  a  slightly  different  form  that  eliminates  the  abso¬ 
lute  value  and  assumes  uniform  wind  and  turbulence.  The  desired  form  is 


Ci(P,vt)  =  <I>*  exp 


-P* 


32  f  V  1^^*  r  f 

~  3  [p  I  {  -  2(1 -w)  PVt  cose  ■+  (Vt)*J  dw 


2ac 


(2) 


where  6  is  the  angle  between  P  and  V.  Differentiating  Eq.(2)  with  respect  to  t  and  set¬ 
ting  the  result  equal  to  zero  yields 


f  ((1— w)^  —  2(1— w)  X  cosH  -P  x^]  1(1— ■"■)  —  x]dw  =  0  (-3) 

0 


•where 


X  — 


(4) 


For  the  case  ■where  V  and  P  are  horizontal  and  V  is  uniform  along  the  path,  Eq.(3)  can 
be  solved  analytically  with  the  result  that  x  =  1/2  or 


V  = 


(5) 


For  the  case  where  P  is  orthogonal  to  V,  the  solution  to  Eq.(3)  is  x  =  0  with  the 
result  that  there  is  no  shift  of  the  peak  of  the  time  delayed  covariance  function  from 
zero  and  consequently  no  sensitivit}'  to  that  component  of  the  wind.  The  relative  wind 
sensitivity  function  for  the  time  delay  to  peak  method  can  be  found  using  Eq.(l)  and 
reformulating  Eq.(3).  Numerical  work  indicated  that  the  relative  sensitivity  function  on 
X  is 


x(w)  =  (1— w)cos6  (6) 

Consequently  the  relative  weighting  on  Tp  the  measured  quantity  is  (1-w),  so  it  would 
appear  that  the  relative  weighting  on  the  measured  wind  along  P  is  l/(l-w).  It  should 
be  noted  that  the  relative  wind  sensitivity  function  is  not  the  conventional  wind  weight¬ 
ing  function.  Because  Eq.(3)  is  a  non-linear  function  of  the  wind,  a  linear  wind  weight¬ 
ing  function  cannot  be  defined.  However,  the  relative  wind  sensitivity  function  does 
give  the  relative  contribution  to  the  result  of  each  part  on  the  path  for  equal  values  of 
wind. 

The  slope  of  the  TLC  function  can  also  be  used  to  process  the  data.  Using  Eq.(  I ) 
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S,  =  -^  C,(P,Vt) 


0 


=  Ci(P,0)  j(5/3)(0.546  Lk^)  f  dw  C^(y,-)  Vp(w)(l-w)2/3  (7) 

•where  Vp  is  the  component  of  -wind  along  P.  From  Eq.(7)  it  can  be  seen  directly  that 
the  wind  weighting  function  is  (1— and  that  the  slope  method  of  processing  is  sen¬ 
sitive  only  to  the  wind  along  P.  Consequently,  ■^’ith  two  orthogonal  sets  of  detectors, 
the  vector  wind  in  the  plane  perpendicular  to  the  line  of  sight  can  be  found.  When  the 
wind  and  turbulence  are  uniform,  Eq.{7)  reduces  to 

qo  p2/3 

Si  =  C,(P,0)  X  “w  ''p 

Po 


which  can  be  used  to  find  Vp  from  measurements  of  the  slope  and  covariance,  Ci(P,0). 
Unfortunately  Eq.(8)  also  depends  on  p^,  which  is  not  known.  Consequently,  the  slope 
method  of  processing  must  be  used  in  conjunction  with  one  of  the  other  methods  in 
order  to  be  useful. 

The  Briggs  method  measures  the  time  delay  at  which  the  autocovariance  and  time 
delayed  covariance  curves  cross.  For  speckle-turbulence  interaction,  a  slightly  modified 
definition  as  illustrated  in  Eq.(9)  is  used. 

exp[PV2ao^)  Ci(P,t)  =  Ci(0,t)  (9) 


Using  Eq.{2)  in  Eq.(9),  the  expression  becomes 

1 

f  1(1— w^)y^  —  2(1— w)  ycos0  -f  1]®/*  dw 
0 


where 


(10) 


y  = 


p 

Vtb 


(11) 


Numerical  work  indicates  that  using 
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y  =  3.1056  cose  (12) 

yields  values  of  y  accurate  to  within  0.05.  Consequently,  the  wind  when  it  is  aligned 

with  P  is  given  by 


3.1056  Tb 

For  the  vector  wind  case,  three  or  more  detectors  are  used  to  form  an  array  where 
there  are  two  vector  spacings  Pj  and  Pj  that  are  orthogonal  (another  possibility  is  three 
detectors  in  an  equilateral  triangle).  In  this  case 


yi  =  - -  =  3.1056  cosfi, 

‘  Vtb,  ' 


y2  =  - -  =  3.1056  cosft2  =  3.1056  sinBj 

Vtb2 


where  Bj  =  6,  —  90° 

Taking  the  ratio  of  Eqs.(14)  and  (15)  yields 


IL  _ 

y2 


=  cotanOj 


From  the  time  delay  measurements  and  the  known  detector  spacings,  the  angle  Bj  can 
be  determined.  The  two  wind  components  are  then  given  by 


VP, 

V,  =  — ^  =  V  cosB,  = 


3.1056  Tb, 


VP2 

V2  =  - -  =  V  sinB,  = 


3.1056  Tb2 


The  relative  wind  sensitivity  function  for  the  Briggs  method  was  not  determined. 
However,  numerical  work  using  a  60%  sinusoidally  modulated  wind  showed  very  lit(!< 
difference  from  the  case  where  the  wind  was  constant. 


mu 
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Another  method  that  may  be  useful  in  processing  the  data  is  the  autocovariance 
width.  In  this  method,  the  time  delay  at  which  the  autocovariance  falls  to  a  prescribed 
value  is  used.  For  speckle-turbulence  interaction,  the  slope  of  the  autocovariance  curve 
is  a  maximum  when 


Ci(o,Vt^)  =  e 


-0.4 


(20) 


which  yields 


V/Po  = 


0.139 


(21) 


As  can  be  seen  from  Eq.(21),  the  strength  of  turbulence  must  be  known  in  order  to  find 
the  wind  using  the  width  of  the  autocovariance  processing  method.  However,  it  can  be 
used  with  other  methods  such  as  the  slope  method  to  find  both  V  and  p^.  In  addition,  it 
should  be  noted  that  this  method  has  no  sensitivity  to  the  direction  of  the  wind.  The 
relative  sensitivity  of  the  width  of  the  autocovariance  method  to  wind  at  different  points 
along  the  path  has  been  determined  numerically  to  be  uniform.  It  is  however,  non¬ 
linear  with  respect  to  the  magnitude  of  the  wind. 

All  of  the  above  methods  have  been  used  to  process  data  collected  by  the  experi¬ 
mental  apparatus.  However,  the  results  obtained  were  not  very  good.  It  is  believed  that 
the  reason  for  this  is  the  rather  small  effect  that  the  turbulence  has  on  the  time  delayed 
covariance  function  at  10.6  microns  wavelength.  As  an  example,  for  a  1000  meter  path- 
length,  10.6  microns  wavelength  and  =  10~*®  the  peak  of  the  TDC  rises  only 

2.5%  above  the  value  at  zero  time  delay.  These  small  changes  coupled  with  having  to 
find  them  in  a  fully  developed  speckle  pattern  makes  processing  difficult.  Howr'^or,  it 
should  be  noted  that  the  above  techniques  should  be  usable  at  much  shorter 
wavelengths  where  the  speckle-turbulence  interaction  is  stronger  and  past  work  in  the 
near  infrared  [4)  tends  to  confirm  that  conclusion. 


Because  of  these  difficulties  with  the  above  techniques  an  effort,  which  turned  out 
to  be  quite  successful,  was  made  to  develop  a  new  processing  technique.  This  technique 
will  be  referred  to  as  the  Z  Log  Ratio  (ZLR)  method.  Two  parameters  B  and  C  are 
defined  as 

B  =  =  (22) 

3  [  Po  J 


^  Ci(P,o)  p2 

C  =  - 1 —  =  exp  -  — 


32  f  P  3 


3  IPo 


where  use  was  made  of  Eq.(2).  Now  using  Eqs.(22)  and  (23)  the  ZLR  is  defined  as 


InC  -I-  F^/2a'‘ 


3 


from  which  the  wind  velocity  can  be  derived  as 


V  =  - 

8 


Z3/6p 


As  will  be  seen  in  a  later  section,  good  results  have  been  achieved  using  the  ZLR 
method.  The  quantities  B  and  C  are  derived  from  the  received  signals,  used  in  Eq.(25) 
to  calculate  Z  which  is  used  in  (26)  to  obtain  the  magnitude  of  the  wind.  It  should  be 
noted  that  there  are  no  arbitrary  parameters  used  in  the  formulation;  and  the  only  sys¬ 
tem  parameters  needed  are  the  spacing  between  the  detectors  and  the  transmitted  beam 
radius.  Since  Eq.(26)  does  not  yield  the  direction  of  the  wind,  for  the  rase  wh'  vr  tlu- 
wind  is  aligned  along  P,  the  sign  of  the  measured  slope  of  the  TLC  can  be  used  for  that 
purpose.  Numerical  analysis  indicates  that  the  wind  sensitivity  function  for  the  ZLR 
processing  technique  is  uniform. 
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The  ZLR  method  may  also  be  used  for  the  vector  wind  case.  Since  Eq.(26)  depends 
on  the  magnitude  of  the  vector  wind,  it  is  still  valid.  All  that  remains  is  to  find  the 
angle  between  P,  and  the  vector  wind.  If  the  slope  is  measured  using  two  sets  of 
orthogonal  detectors  with  spacings  Pj  and  Pj,  then  the  ratio  of  the  slopes  yields 


S,  Ci(P„o)  fP, 


.2/3 


Cl(P2,o)  IP2 


—  j  cotanBj 


(27) 


which  reduces  to 


Si 

—  =  cotan0, 

$2 

for  Pj  =  P2.  The  two  components  of  the  wind  then  are  given  by 

Vj  =  V  cosHj 

and 


(28) 


(29) 


V2  ==  V  sinGj 

In  practice,  Z  is  measured  at  several  time  delays  and  the  results  averaged, 
done  then 


(30) 
If  this  is 


Ivl 


/ 

8 


3/5  M 
'  p  N 

—  V 

N  .ti 


(31) 


where  is  the  measured  value  of  Z  at  time  delay  Tn-  Most  of  the  wind  data  shown  in 
the  next  section  was  processed  using  Eq.(31). 

In  an  attempt  to  further  improve  the  ZLR  method,  it  has  been  reformulated  into  a 
result  that  minimizes  the  mean  square  error  on  the  measured  quantity  Zj,.  The  mean 
square  error  is  given  by 


Using 


MSE  = 


N 


N 

V 

n=l 


*n 

5/6  ■ 

3 

,  p^  . 

(32) 
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aMSE 


the  wind  velocity  that  minimizes  the  mean  square  error  is  given  by 


N 

3/8  V 

n=l 


^MMSE 


Now  putting  Eq.(34)  back  into  Eq.(32)  the  minimum  mean  square  error  is  found  as 


1  r  N  f  N  N 

MMSE=^  2  z^-  V  /V  ItJ 

N  n=i  D=1  0=1 


The  use  of  Eq.(34)  usually  produces  slightly  better  results  than  Eq.(31).  It  is  proposed 


that  Eq.(35)  might  be  used  as  a  "goodness*  measure  on  the  data. 


In  a  subsequent  section  the  experimental  results  using  Eq.(31)  and  Eq.(34)  will  be 


presented  along  with  one  more  variation  on  the  processing  that  gives  both  a  better  esti¬ 


mate  of  the  wind  and  significantly  reduces  the  statistical  fluctuation  in  the  measure¬ 


ment.  The  data  is  preprocessed  to  yield  a  one  bit  per  sample  binary  data  set.  This  is 


accomplished  by  letting  the  bit  be  one  or  zero  depending  on  whether  the  received  inten¬ 


sity  is  greater  than  or  less  than  the  mean  respectively.  This  also  significantly  decreases 
the  processing  time  required  and  will  make  designing  a  portable  instrument  much 


easier.  This  technique  will  be  referred  to  as  Binary-Z  Log  Ratio  (B-ZLR)  processing. 


Experimental  Design 


The  transmitter/receiver  optical  design  is  shown  in  Figure  3.  The  source  is  n 


waveguide  laser  operating  at  10.6  microns  wavelength.  The  laser  beam  is  first  expanded 


by  a  factor  of  3  and  then  most  of  the  beam  is  split  off  by  a  90'^  beam  splitter  an 
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directed  through  an  acoustooptic  modulator  (AOM)  used  as  a  frequency  shifter  and 
driven  at  37.5  MHz.  The  frequency  shifted  beam  at  a  little  over  one  watt  of  power  then 
passes  through  a  quarter-wave  plate  to  produce  a  circularly  polarized  beam  that  is 
focused  onto  the  target  by  a  10  X  beam  expander.  The  half-wave  plate  at  the  output  of 
the  3  X  beam  expander  is  used  to  change  the  polarization  of  the  laser  beam  from  verti¬ 
cal  to  horizontal  to  match  the  requirements  of  the  AOM.  A  circularly  polarized  beam  is 
transmitted  in  order  to  eliminate  fluctuations  in  the  returned  signal  due  to  target 
induced  polarization  changes. 

Approximately  10%  of  the  laser  beam  passes  through  the  90%  beam  splitter  for 
use  as  an  optical  local  oscillator.  It  passes  through  a  second  AOM  which  is  also  used  as 
a  frequency  shifter  and  driven  at  42.5  MHz.  The  purpose  of  the  second  AOM  is  to  pro¬ 
vide  isolation  between  the  transmitted  beam  and  the  local  oscillator.  An  isolation  of 
over  200  db  has  been  achieved  [5]  which  is  sufficient  not  to  degrade  the  extreme  sensi¬ 
tivity  of  optical  heterodyne  detection.  The  frequency  shifted  local  oscillator  beam  is 
then  circularly  polarized  by  a  quarter-wave  plate  to  match  the  transmitted  beam,  split 
into  two  parts  and  focused  onto  two  HgCdTe  detectors.  The  detectors  are  photovoltaic, 
0.2  mm  square  and  were  housed  in  separate  dewars  to  allow  some  flexibility  in  selecting 
the  spacing  between  the  detectors  and  to  simplify  the  optical  design. 

The  transmitted  beam  is  scattered  diffusely  by  a  4  foot  x  4  foot  sandblasted  alumi¬ 
num  target.  The  speckle  field  that  is  generated  by  the  target  and  modulated  by  the 
turbulence  in  the  atmosphere  returns  to  the  transmitter /receiver  where  it  is  directed  by 
two,  one-inch  mirrors  set  at  45  degrees  onto  lenses  that  focus  it  onto  the  dpterfoi'  TIm 
returned  signal  and  the  optical  local  oscillator  photomix  at  the  detector  to  produce  ;i 
heterodyne  signal  at  their  difference  frequency  of  5  MHz.  This  signal  is  then  processed 
to  obtain  the  wind.  Since  the  detectors  are  spaced  horizontally,  the  system  meacur'-'' 
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the  horizontal  component  of  the  wind.  If  the  sector  wind  were  to  be  measured,  then  a 
minimum  of  one  more  detector,  spaced  vertically,  would  be  required. 

The  receiver  electronics  for  one  channel  is  illustrated  in  Figure  4.  The  heterodyne 
signal  current  at  5  MHz  from  each  detector  is  amplified,  filtered,  mixed  down  to  100 
kHz,  precision  rectified  to  recover  the  field  amplitudes,  and  low  pass  filtered  before 
being  digitized  for  processing.  Because  the  interaction  with  the  turbulence  is  small, 
good  linearity  (0.1%)  and  stability  are  required.  The  received  intensities  are  recovered 
by  squaring  the  measured  amplitudes  in  a  computer. 

Experimental  Results 

Measurements  were  made  using  500  meter  and  1000  meter  paths  over  very  flat, 
featureless  terrain  where  the  turbulence  is  quite  uniform  over  the  path.  In  situ  data 
was  taken  using  a  Campbell  Scientific  CA-9  path  averaging  scintillometer  and  crosswind 
sensor.  The  Campbell  unit  is  a  double  ended  unit  with  a  transmitter  at  one  end  and  a 
receiver  at  the  other  and  has  a  cosine-like  wind  sensitivity  function;  whereas  the  hetero¬ 
dyne  system,  utilizing  the  ZLR  method  of  processing  has  a  uniform  wind  sensitivity 
function.  Since  the  wind  was  not  uniform  along  the  path,  some  of  the  differences 
between  the  heterodyne  measurements  and  the  in  situ  data  may  be  attributed  to  the 
difference  in  their  wind  sensitivity  functions.  Figures  5  and  6  show  wind  magnitude 
data  taken  with  the  target  at  500  and  1000  meters,  with  approximately  2  and  4  meters 
per  second  wind  respectively  and  at  an  intermediate  level  of  turbulence.  The  overall 
results  are  good.  There  is  moderate  statistical  fluctuation  about  the  me.an  du' 
speckle  and  short  averaging  times;  but  the  mean  values  averaged  over  the  100  s‘'r'nnds 
of  data  are  within  0.12  and  0.54  m/sec  respectively  of  the  in  situ  measurements.  It 
should  be  noted  that  the  fluctuation  can  be  reduced  either  by  longer  time  averac<‘' 


and/or  by  using  additional  detectors  and  spatial  averaging.  In  addition,  as  will  be 
demostrated  later  in  his  section,  the  use  of  B-ZLR  or  MMSE  processing  also  reduces  the 
fluctuation  due  to  speckle.  The  wind  direction  was  determined  by  the  sign  of  the  slope, 
for  the  second  averages,  the  correct  sign  of  the  slope  was  measured  about  70%  of  the 
time.  This  result  is  not  as  good  as  was  expected;  and  work  is  continuing  to  improve  the 
processing. 

An  important  consideration  in  using  speckle-turbulence  interaction  for  measuring 
winds  is  operation  at  low  turbulence  levels.  An  estimate  of  the  lower  limit  at  which  this 
type  of  system  should  operate  will  be  made,  followed  by  some  experimental  work  at  low 
turbulence. 

Aji  estimate  of  the  sensitivity  of  the  experimental  system  can  be  made  by  consider¬ 
ing  Eq.(8).  Since  the  slope  is  proportional  to  Cq.  as  decreases  the  slope  decreases 
and  eventually  a  limit  will  be  reached  below  which  the  system  will  not  respond.  Also 
there  are  deterministic  (as  opposed  to  random)  errors  such  as  gain  drifts  for  the  two 
optical  channels  that  do  not  average  to  zero.  Consequently,  it  is  doubtful  that  the  sys¬ 
tem  can  measure  the  received  intensity  normalized  by  the  mean  to  better  than  one- 
tenth  of  one  percent.  This  translates  to  no  worse  than  four-tenths  of  one  percent  on  the 
normalized  slope.  Using  this  with  Eq.(8),  it  is  found  that 

2  >  1-741  X 
“  ~  Lp2/3VpC,(P) 

It  can  be  seen  that  as  the  wavelength  \  decreases  and  the  path  length  increases,  the  sen¬ 
sitivity  improves.  For  the  experimental  system  under  consideration  and  with  L  —  K'nn 
m  and  Vp  =  1  m/s,  Eq.(36)  becomes 


This  is  the  minimum  C„  at  which  the  system  is  expected  to  respond  under  the  condi¬ 
tions  given  and  appears  to  be  low  enough  to  make  the  proposed  technique  quite  useful. 
It  is  expected  then,  that  a  system  of  this  type  should  provide  good  measurements  at 
levels  around  10  times  the  level  predicted  by  (37)  or  3  x  10“*® 

Figure  7  shows  a  data  set  processed  using  the  ZLR  method  and  corresponding  to 
C„  =  7.4  X  10“*®  m“^^  which  constitutes  a  low  value  of  turbulence.  There  is  the  usual 
fluctuation  about  the  mean  and  the  mean  value  of  1.64  m/sec  averaged  over  the  100 
seconds  of  the  data  set,  is  less  than  the  mean  value  of  1.96  m/sec  measured  by  the 
Campbell  unit.  This  reduction  in  mean  value  is  believed  to  be  indicative  that  the  system 
is  just  beginning  to  lose  sensitivity  (with  ZLR  processing). 

Figure  8  shows  this  same  data  set  processed  using  the  MMSE  method  and  it  can  be 
seen  that  the  error  in  the  mean  is  comparable  to  the  ZLR  method;  but  the  statistical 
fluctuations  are  reduced  yielding  an  rms  error  of  0.56  m/sec  as  opposed  to  0.79  for  ZLR 
processing. 

The  results  of  processing  this  data  using  the  B-ZLR  method  are  shown  in  Figure  9. 
There  is  a  substantial  reduction  in  error  and  fluctuation  as  compared  to  the  results 
obtained  by  ZLR  processing.  The  mean  error  was  only  0.15  m/sec  and  the  RMS  error 
only  0.35  m/sec.  Processing  the  data  using  the  binary  data  set  and  the  MMSE  method 
did  not  yield  results  any  better  than  the  B-ZLR  method.  Considering  that  the  MMSE 
method  requires  more  processing  than  the  ZLR  method  it  is  concluded  that  the  best  of 
the  methods  that  were  considered  is  the  B-ZLR  method  of  processing  and  that  good 
measurements  can  be  made  at  turbulence  levels  down  into  the  10“’®  ^  rang* 
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Figure  2  -  Processing  Techniques 


RMS  Error  =  0.56  m/s 


